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Abstract. Wc investigate polarization of high-energy 
emissions from the Crab pulsar in the frame work of the 
outer gap accelerator. The recent version of the outer gap, 
which extends from inside the null charge surface to the 
light cylinder, is used for examining the light curve, the 
spectrum and the polarization characteristics, simultane- 
ously The polarization position angle curve and the polar- 
ization degree are calculated to compare with the Crab op- 
tical data. We show that the outer gap model explains the 
general features of the observed light curve, the spectrum 
and the polarization by taking into account the emissions 
from inside of the null charge surface and from tertiary 
pairs, which were produced by the high-energy photons 
from the secondary pairs. For the Crab pulsar, the po- 
larization position angle curve indicates that the viewing 
angle of the observer measured from the rotational axis is 
greater than 90°. 



1. Introduction 

Young pulsars such as the Crab pulsar are strong 7-ray 
sources. The EGRET instrument revealed that the light 
curve with double peaks in a period and the spectrum 
extending to above GeV are typical features of the high- 
energy emissions from the 7-ray pulsars. Although these 
data have constrained proposed models, the origin of the 
7-ray emission is not yet conclusive. On important reason 
is that various models have successfully explained the fea- 
tures of the observed spectra and/or light curves. For ex- 
ample, the polar cap model (Daugherty & Harding 1996), 
the caustic model (Dyks et al. 2004) and the outer gap 
model (Cheng et al. 2000, hereafter CRZ00), all expect 
the main features of the observed light curve. So, we can- 
not discriminate the three different models using the light 
curve. Furthermore, both polar cap and outer gap models 
have explained the observed 7-ray spectrum (Daugherty 
& Harding 1996; Romani 1996). 

Polarization measurement will play an important role 
to discriminate the various models, because it increases 
observed parameters, namely, polarization degree (p.d.) 
and position angle (p. a.) swing. So far, only the optical 



polarization data for the Crab pulsar is available (Kan- 
bach et al. 2005) in high energy bands. For the Crab pul- 
sar, the spectrum is continuously extending from optical 
to 7-ray bands. In addition, their pulse positions are con- 
sidered so well that the optical emission mechanism is re- 
lated to higher energy emission mechanisms. In the future, 
the next generation Compton telescope will probably be 
able to measure polarization characteristics in MeV bands. 
These data will be useful for discriminating the different 
models. 

In this paper, we examine the optical polarization char- 
acteristics of the Crab pulsar with the light curve and the 
spectrum in frame works of the outer gap model. CRZ00 
has calculated the synchrotron self-inverse Compton scat- 
tering process of the secondary pairs produced outside the 
outer gap and has explained the Crab spectrum from X- 
ray to 7-ray bands. In CRZ00, however, the outer-wing 
and the off-pulse emissions of the Crab pulsars cannot 
be reproduced, because the traditional outer gap geom- 
etry, which extends from the null charge surface of the 
Goldrcich-Julian charge density to the light cylinder, is 
assumed. Furthermore, the spectrum in the optical band 
was not considered. In this paper, on these grounds, we 
modify the CRZ00 geometrical model into a more realis- 
tic model, following recent 2-D electrodynamical studies 
(Takata et al. 2004; Hirotani 2006), and calculate the light 
curve, the spectrum and the polarization characteristics of 
the Crab pulsar. 

2. Calculation method 

The outline of the outer gap model for the Crab pul- 
sar is as follows. The charge particles are accelerated by 
the electric field parallel to the magnetic field lines in so 
called gap, where the charge density is different from the 
Goldreich-Julina charge density. The high energy particles 
accelerated in the gap emit the 7-ray photons (called pri- 
mary photons) via the curvature radiation process. For the 
Crab pulsar, most of the primary photons escaping from 
the outer gap will convert into secondary pairs outside the 
gap, where the accelerating electric field vanishes, by col- 
liding with synchrotron X-rays emitted by the secondary 
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pairs. The secondary pairs emit optical - MeV photons via 
the synchrotron process and photons above MeV with the 
inverse Compton process. The high-energy photons emit- 
ted by the secondary pairs may convert into tertiary pairs 
at higher altitude by colliding with the soft X-rays from 
the stellar surface. The tertiary pairs emit the optical-UV 
photons via the synchrotron process. This secondary and 
tertiary photons appear as the observed radiations from 
the Crab pulsar 

2.1. Outer gap 

We consider the outer gap geometry that is extending from 
inside null charge surface to near the light cylinder (Takata 
et al 2004; Hirotani 2006). Because the Crab pulsar has a 
thin gap, we describe the accelerating electric field (Cheng 
ct al. 1986) with 



Eu(r) 



nB(r)f 2 (r)Rj 
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where f(r) is the local gap thickness in units of the light 
radius, Ri c = c/f2, and s(r) is the curvature radius of the 
magnetic field line. The typical fractional size of the outer 
gap is given by f{Ri c /2) ~ 5.5B^ 4 / 7 P 26 / 21 (~ 0.04 for 
the Crab pulsar) and the local fractional size is estimated 
by f(r) ~ /(i? ic /2)(2r/i? ic ) 1 - 5 (Zhang & Cheng 1997). 

The electric field described by equation ([1]) can be 
adopted for the acceleration beyond the null charge sur- 
face. Inside null charge surface, the electric field rapidly 
decreases because of the screening effects of the pairs. To 
simulate the accelerating electric field inside null charge 
surface, we assume the following field, 
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n < r < r n , 
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where E n is the strength of the electric field at the null 
charge surface and r n and are the radial distances to 
the null charge surface and the inner boundary of the 
gap, respectively. The local Lorcntz factor of the accel- 
erated particles in the outer gap is described by T e (r) = 
[3s 2 (r)£'j|/2e] 1 / 4 with assuming the force balance between 
the acceleration and the curvature radiation back reaction. 

We assume that the outer gap extends around the 
whole polar cap. In the calculation, we constrain the 
boundaries of the axial distance and radial distance for 
the emission regions with p max = 0.9Ri c and r = Ri c , 
respectively. 

2.2. Synchrotron radiation from the pairs 

We calculate the synchrotron radiation at each radiating 
point following CRZ00. The photon spectrum of the syn- 



chrotron radiation is (CRZ00) 

3 1 / 2 e 3 B(r) sin p (r) 



Esyn (,Ej , r) 



mc 2 hE~, 



dn e (r) 



F{x)dE e dV rad , 



(3) 



where x = Ey/E syn , E syn is the typical photon energy, 9 p 
is the pitch angle of the particle, F(x) = x K 5 / 3 (y)dy, 
where K 5 / 3 is the modified Bessel function of order 5/3, 
dn e /dE e cx E~ 2 is the distribution of the pairs and dV ra d 
is the volume element of the radiation region considered. 

The pitch angle of the secondary pairs is estimated 
from sin# p (i?z c ) ~ A/s(i?; c ), where A is the mean free 
path of the pair-creation between the primary 7-rays 
and the non-thermal X-rays from the secondary pairs. 
For the Crab pulsar, the mean free path becomes A ~ 
(rix<T 71 ) _1 ~ 10 7 cm, where we used the typical num- 
ber density n x ~ L x (< E x >)/5QRf c c < E x >~ 
8 x 10 17 cm~ 3 , the luminosity L x ~ 10 35 erg/s, the typ- 
ical soft-photon energy for the pair-creation, < E x >^ 
(2m e c 2 ) 2 /10 GeV ~ 100 eV, and cr 77 ~ cr T /3. There- 
fore the pitch angle of the secondary pairs at the light 
cylinder is estimated as sin^p ~ X/s(Ri c ) ~ 0.06, 
and the local pitch angle is calculated from sin 6> p (r) = 
sm6 p {R lc ){r/R lc ) 1 ^. 

Some high-energy photons emitted by the inverse 
Compton process of the secondary pairs may convert into 
tertiary pairs at higher altitude by colliding with thermal 
X-ray photons from the star. In this paper, we assume 
that the maximum energy of and the local number den- 
sity of the tertiary pairs are smaller than about 10% of 
those of the secondary pairs. Because the pitch angle of 
the pairs increases with altitude, we use sin# p = 0.1 for 
the pitch angle of the tertiary pairs. In fact, the results 
are not sensitive to the pitch angle of the tertiary pairs. 



2.3. Stokes parameters 

For a high Lorentz factor, we can anticipate that the emis- 
sion direction of the particles coincides with the direction 
of the particle's velocity. In the inertial observer frame, 
the particle motion may be described by 



A> 

cos 9pb + (3q sin 9 p b± + /3 co e0, 



(4) 



where the first term in the right hand side represents the 
particle motion along the field line, for which we use the 
rotating dipole field, and b is the unit vector of the mag- 
netic field line. The second term in equation ([4]) repre- 
sents gyration motion around the magnetic field line and 
b± = cos Scj)k + sin 8(j>k x b is the unit vector perpendicu- 
lar to the magnetic field line, where 54> refers the phase of 
gyration motion and k = (6- V)6/|(b- V)6| is the unit vec- 
tor of the curvature of the magnetic field lines. The third 
term is co-rotation motion with the star, f3 co = p^l/c. The 
emission direction of equation (j4|) is described in terms 
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of the viewing angle measured from the rotational axis, 
£ = cos -1 n z , and the rotation phase, $ = — $„ — r ■ n, 
where n z is the component of the emission direction par- 
allel to the rotational axis, 4>„ is the azimuthal angle of 
the emission direction and r is the emitting location in 
units of the light radius. 

Because the particles distribute on the gyration phase 
8(f), the emitted beam at each point must become cone 
like shape with opening angle p {r). We calculate the ra- 
diations from the particles for all of the gyration phase 
5(f) = 2ni/n (i = 1, ■ • ■ , n — 1). 

We assume that the radiation at each point linearly 
polarizes with degree of U syn = (p+l)/(p + 7/3), where 
p is the power law index of the particle distribution, and 
circular polarization is zero, that is, V = in terms of 
the Stokes parameters. The direction of the electric vec- 
tor of the electro-magnetic wave toward the observer is 
parallel to the projected direction of the acceleration of 
the particle on the sky, that is, E em oc a — (n ■ a)n 
(Blaskicwicz et al. 1991). In the present case, the ac- 
celeration with equation (j4|) is approximately written by 
a ~ Pqlob sin# p (— smS(f>k + cos 8<pk x b). 

We define the position angle % 4 to be angle between 
the electric field E em and the projected rotational axis 
on the sky. The Stokes parameters Q l and U l at each ra- 
diation is represented by Q % = H syn P cos 2\ l and U l = 
H syn P sin 2x l - After collecting the photons from the pos- 
sible points for each rotation phase <f> and a viewing angle 
£, the expected p.d. and p. a. are, respectively, obtained 
from P(£,$) = n s „ n> /Q2(£,$)W 2 (£, and 
X (£,$) =0.5atan [[/(£, $)/Q(£,$)]. 

The inclination angle a and the viewing angles £ mea- 
sured from the rotational axis are the model parameters. 
The radial distance of the inner boundary in equation 
(0) is also a model parameter, because the distance rj is 
determined by the current through the gap (Takata et al. 
2004). Because the last-open line must be modified from 
the traditional magnetic surface, which is tangent to the 
light cylinder for the vacuum case, by the plasma effects 
(Romani 1996), the altitude of the upper surface of the 
outer gap, above which the pairs are created and emit the 
synchrotron photons, is used as a model parameter using 
a fractional polar angle a = 9 u /9i c , where 9 U is the po- 
lar angle of the footpoints of the magnetic field lines of 
the gap upper surface and 9i c is the polar angle of the 
field lines which are tangent to the light cylinder for the 
vacuum case. 

3. Results 

Figure [T] compares the polarization characteristics at 1 eV 
predicted by the traditional (left column) and the present 
(right column) outer gap models. The traditional model 
considers the emissions from the secondary pairs with the 
outer gap starting from the null charge surface of the 
Goldrcich- Julian charge density, = r n . The present 
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Fig. 1. Polarization characteristics for the traditional 
model (left column) and the present model (right column) . 
The upper, middle and lower panels in each column show, 
respectively, the light curve, the position angle and the 
polarization degree for a = 40° and a = 0.93. 



model takes into account the emissions from inside null 
charge surface and the tertiary pairs. We assume the outer 
gap starts from the radial distance of 68% of the distance 
to null charge surface, r, = 0.68r„. The model parameters 
are a = 40°, a = 0.93 and £ ~ 101°, where the viewing 
angle is chosen so that the predicted phase separation be- 
tween the two peaks is consistent with the observed value 
<5$ ~ 0.4 phase. 

From the pulse profiles, we find that the traditional 
model can not produce the outer-wing and the off-pulse 
emissions. On the other hand, the present model produces 
the outer-wing and the off-pulse emissions with the emis- 
sions from inside of the null charge surface. 

As seen in the polarization characteristics by the tradi- 
tional model, we find that the secondary emissions beyond 
the null charge surface make the polarization characteris- 
tics such that the polarization degree takes a lower value at 
the bridge phase and a larger value near the peaks. In the 
synchrotron case, the cone like beam is radiated at each 
point, and an overlap of the radiations from the different 
particles on the gyration phase causes the depolarization. 
For the viewing angle £ ~ 101°, the radiations from all 
gyration phases contribute to the observed radiation at 
the bridge phase, because the line of sight passes through 
middle part of the emission regions. In such a case, the 
depolarization is strong, and as a result, the emerging ra- 
diation from the secondary pairs polarizes with a very low 
p.d. (< 10%). Near the peaks, the radiations from the 
some range of gyration phases are not observed because 
the line of sight passes through the edge of the emission re- 
gions at peaks. In such a case, the depolarization is weaker 
and the emerging radiation highly polarizes. 

In the polarization angle swing in the present model, 
the large swings appear at the both peaks, and the differ- 
ence of the position angle between the off-pulse and the 
bridge phases is about 90 degree. We can see the effects 
of the tertiary pairs on the polarization characteristics at 
the bridge phase. By comparing the p.d. at the bridge 
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A. 




Fig. 2. Optical polarization for the Crab pulsar. 
Left:Polarization characteristics for the total emissions. 
Right :Polarization characteristics of the emissions after 
subtraction of the DC level (Kanbach et al. 2005). The 
figures was transcribed from Dyks ct al. (2004). 



Fig. 3. The optical-X ray spectrum for the Crab pulsar. 
The calculation is for a = 40°, a = 0.93, r\ = 0.68 and 
£ ~ 101°. The X-ray data are taken from Kuiper et al. 
(2002) and reference therein, and the optical data from 
Sollerman et al. (2000). 



phase, we find that tertiary pairs produce the radiations 
with ~ 10% of the p.d. at the bridge phase. 

Figure [2] summarizes the Crab optical data. Left col- 
umn and the right column show, respectively, the Crab 
optical data for the total emissions and for the emissions 
after subtraction of the DC level, which has the constant 
intensity at the level of 1.24% of the main pulse intensity. 
In the total emissions (left column) , the impressive polar- 
ization feature is that the off-pulse and bridge phases have 
the fixed value of the p. a. These polarization features of 
the observation are not predicted by the present model. 
The present model are more consistent with the Crab op- 
tical data after the subtraction of the DC level. Especially, 
the model reproduces the most striking feature in the ob- 
served p. a. that the large swing at both peaks, and the 
observed low p.d. at bridge phase ~ 10%. Also, the pat- 
tern of the p.d. are reproduced by the present model. 

Figure [3] compares the model spectrum with the Crab 
data in optical-MeV bands. The model parameters are 
same with that in the right column in Figure [2] In this 
case, we assume that the pairs escape from the light cylin- 
der with the Lorentz factor r ~ 20 due to the synchrotron 
cooling effect, which makes a spectral break around 10 eV 
in Figure [3] The model spectrum also explains the general 
features of the data. The outer gap model can explain the 
general features of the observed light curve, the spectrum 
and the polarization characteristics in optical band for the 
Crab pulsar, simultaneously. 

It is worth to note that we can distinguish the two 
viewing angles mutually symmetric with respect to the 
rotational equator using the polarization swing. For such 
symmetric viewing angles, the light curves, the spectra 
and the p.d. curves are identical. However, the p. a. curves 
arc mirror symmetry with respect to the rotational equa- 
tor because of the difference directions of the projected 
magnetic field on the sky. The pattern of the position an- 
gle for the viewing angle smaller than 90° swings to op- 
posite directions from the Crab data at the both peaks. 



Therefore, the present model predicts that the viewing 
angle larger than 90° are preferred for the Crab pulsar. 

4. Conclusions 

Wc have considered the light curve, the spectrum and the 
polarization characteristics for the Crab pulsar predicted 
by the outer gap model which takes into account the emis- 
sions from inside the null charge surface and from the ter- 
tiary pairs. We have shown that the expected polarization 
characteristics are consistent with the Crab optical data 
after subtraction of the DC level. The outer gap model 
explains the spectrum, light curve and the polarization 
characteristics , simultaneously. 
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